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ABSTRACT 

We study the properties of a significant thermal emission component that was identified in 47 GRBs 
observed by different instruments. Within the framework of the “fireball” model, we deduce the values 
of the Lorentz factor P, and the acceleration radius, r o, for these bursts. We find that all the values 
of r in our sample are in the range 10 2 < T < 10 3 , with (r) = 370. We find a very weak dependence 
of P on the acceleration radius ro, T oc rff with a = —0.10 ± 0.09 at the a = 2.1 confidence level. 

The values of ro span a wide range, 10 6 ' 5 < ro < 10 9 5 cm, with a mean value of (ro) ~ 10 s cm. This 
is higher than the gravitational radius of a 10 Mq black hole by a factor « 30. We argue that this 
result provides indirect evidence for jet propagation inside a massive star, and suggests the existence 
of recollimation shocks that take place close to this radius. 


1. INTRODUCTION 

One of the major developments in the study of gamma- 
ray bursts (GRBs) in recent years has been the realiza¬ 
tion that a thermal component may be a key spectral 
ingredient. While the shape of most GRB spectra do not 
resemble a “Planck” function, in a non-negligible minor¬ 
ity of GRBs a careful spectral analysis reveals a spec¬ 
tral component that is consistent with having a black- 
body (“Planck”) shape, accompanied by an additional, 
non-t hermal part. Following pioneering work by iRvdel 
(|200412005H . such a componen t was clearly identifi ed in 
56 GRBs detected by BATSE (|Rvde fc Pe’erl[2009tl and 
is now identified in several Fermi GRB s as we ll. A few 
notable examp les are GRB090902B (iR.vde et all [201fll l. 
GRB100507 (IGhirla/nda et a.1.1 1201 .4 . GRB100724B 

(IGuiriec et all 1201 111. GRB110721A (lAxelsson et al. 

120121 llvvani et all 1201311 . GRB120323A (IGuiriec et,~ 
I20TT and GRB101219B (ILarsson et a1J[20T5h FI 

The existence of a thermal emission component may 
hold the key to understanding the prompt emission spec¬ 
tra. First, this component provides a physical explana¬ 
tion to at least part of the observed spectra. Further¬ 
more, thermal photons serve as seed photons for inverse 
Compton (IC) scattering by energetic electrons, therefore 
these photons may play an important role in explaining 
the non-thermal part of the spectra as well. In f act, a s 
it was recentl y shown by lAxelsson &; Borgonovol (|2015l) ; 
lYu et all (|2015bf h while clear “Planck” spectra are only 
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rarely observed, the narrowness of the spectral width of 
GRBs rules out a pure synchrotron origin in nearly 100% 
of the GRB spectra observed to date. Thus, it is pos¬ 
sible that modified “Planck” spectra contribute to the 
observed emission in a very large fraction of GRBs. 

Thermal photons decouple from the plasma at the pho¬ 
tosphere, which is by definition the inner most region 
from which an electromagnetic signal can reach the ob¬ 
server. Therefore the properties of a Plank spectral com¬ 
ponent directly reveal the physical conditions at the pho¬ 
tosphere, in those GRBs in which it can be directly iden¬ 
tified. This feature is in contrast to the non-thermal 
spectral component, the origin of which is still uncer¬ 
tain. This is due to the fact that the radiative origin 
of the non-thermal component is still debatable, and its 
exact emission radius is very poorly constrained, both 
theoretically and observationally. 

In the framework of the classical “fireba ll” model 
(|Paczvnskilll986t iRees &; Meszarodl 1 9921 179941) . the pho- 
tospheric radius, r p h depends only on two free model 
parameters: the luminosity, L , and the L orentz fac¬ 
tor, r at the photospheric r adiujU (e.g.. iPaczvnskil 
1990t lAbramowicz et al.l 119911 IMeszdros fe Reesl 120001 
MeszaroslI2006L and references therein). The observed 
temperature weakly depends, in addition, on a third pa¬ 
rameter, ro, as T ob oc rj/ 6 . Here, ro is the acceleration ra¬ 
dius, which is the radius where the acceleration of plasma 
to relativistic (kinetic) motion begins. Thus, by defini¬ 
tion, at this radius the bulk Lorentz factor T(ro) = 1, 
while at larger radii T(r) increases at the expense of the 
internal energy. In the classical “fireball” model, where 
the outflow expands freely and magnetic fields are sub¬ 
dominant, the growth is linear, T(r) oc r below the sat¬ 
uration radius r sat , and T(r) oc r° above this radius, as 
all the available energy is already in kinetic form. 

The simple dependences of the temperature and photo- 
spheric radius on the luminosity, T and rg, have a strong 

7 This statement holds under the assumption that the photo- 
spheric radius r p h is larger than the saturation radius, r sa t, which 
is the radius at which all the available internal energy is converted 
to kinetic energy. 
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implication. As was shown bv iPe’er et al.l (120071) . one 
can use the measured values of the temperature and ob¬ 
served flux to directly measure r p h and T. Using the 
fireball model scaling laws, one can then infer the value 
of ?'o- Therefore, using only observational quantities one 
can deduce the entire (basic) dynamics of the earliest 
stages of the jet evolution within the framework of the 
“fireball” model. These dynamics, in turn, provide a 
strong tool in constraining models of GRB progenitors. 

In order to perform these calculations, knowledge of 
the distance to the GRB is required, as the observed flux 
needs to be converted to luminosity. Moreover, reliable 
estimates of the dynamical parameters (r, r p h and ro) 
rely on the assumption that the observed thermal compo¬ 
nent is not strongl y distorted (e.g ., by sub-photosplieric 
dissipation; see IPe’er et al.l 120061 ). Validating this last 
assumption is, nonetheless, relatively easy, as a signifi¬ 
cant thermal component F t h 1$ F 1 that can be directly 
observed, necessitates tha t a stro ng distortion does not 
occur (see, however. IAhlgren et alJl2()T5l) . 

Unfortunately, as of now, there are only very few GRBs 
which fulfill both requirements, namely, (1) a significant 
thermal component is clearly observed in their spectra, 
and (2) their luminosity distance is known. These obser¬ 
vational constraints limited, so far, the ability to carry 
out a statistical study of the outflow properties derived 
from photospheric emission. 

Nevertheless, in recent years the number of GRBs in 
which a thermal component could be clearly identified is 
rapidly increasing, and is now at a few dozens. Further¬ 
more, as we point out in the present work, uncertainties 
in the redshift have, in fact, only a weak effect on the 
deduced value of ro. Consequentially, a good estimate of 
the value of ro can be achieved even for those GRBs for 
which the redshift is unknown. This enables the study of 
the hydrodynamical parameters in a reasonable sample 
of GRBs for which a thermal component was identified. 

In this paper we analyze the existing data for those 
GRBs in which a distinct thermal component was clearly 
identified. Using the existing data, we deduce the values 
of the bulk Lorentz factor, T, and the acceleration radius, 
ro- We find an average value of T ~ 10 25 , which is simi¬ 
lar to previous estimates, based on various methods. On 
the other hand, we find that (ro) ~ 10 8 cm, higher than 
the common assumption of < 10' cm. These larger than 
expected ro values may be due to propagation effects of 
the jet, such as recollimation shocks, mass entrainment, 
or a non conical structure, and may provide indirect ev¬ 
idence for the presence of a massive progenitor star. 

This paper is organized as follows. In section [2] we 
describe our sample selection and method of analysis. 
Our results are presented in section [3] We discuss our 
findings in section [51 before summarizing and concluding 
in section [5] 

2. SAMPLE AND METHOD OF ANALYSIS 


flow within the framework of the “fireball” model 0 via 
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Here, or is Thomson’s cross section, 1Z = 
(Fg B /aT ob4 ) 112 , Fg B is the observed black body 
flux, T ob is the temperature of the thermal component, 
and a is Stefan-Boltzmann constant. The total (thermal 
+ non-thermal) flux is denoted by F ob , and Y > 1 is the 
ratio between the total energy released in the explosion 
producing the GRB and the energy observed in y-rays 
(both thermal and non-thermal). Measurements of 
the exact value of Y are difficult to conduct, though 
estimates ca n be d one using a ftergl o w observations 
(ICenko et al.1 120111 : IPe’er et all 120121 IWveoda et al.1 
1201, ’ill . 

A key difficulty in identifying a thermal component 
and conducting the calculations is the fact that the sig¬ 
nal varies with time. Thus, a time_ depende nt an alysis 
is required. As was shown by IRvdel (|2004 1200511 and 
IRvde fe Pe’eil (|2009H . for bursts which show relatively 
smooth, long pulses, and in which a thermal component 
could be identified, both the temperature and the flux 
show a very typical behavior: a broken power law in 
time. Typically, before t% b k ~ a few s, the tempera¬ 
ture is roughly constant, while the thermal flux increases 
roughly as F B b B oc t 1 / 3 . At later times, both the temper¬ 
ature and the flux decay, T ob oc t -2 / 3 and F B b B oc t~ 2 . 
Although the non-thermal and thermal fluxes are often 
correlated, the break time does not always coincide with 
the peak of the (non-thermal) flux. In many GRBs it is 
associ ated wit h the begin ning of the rapid decay of the 
pulse (IRvde fc Pe’efi[2009i) . 0 

The exact interpretation of the early and late tem¬ 
poral behavior is still not fully clear. An interesting 
finding is that during the rise phase of the pulses (the 
first few seconds), the temperature is roughly constant. 
This implies via Equation [T] that the Lorentz factor is 
roughly constant (assuming that the ratio Fbb/F tot does 
not vary much), though the total flux changes substan¬ 
tially. This no longer holds during the decay phase. A 
leading idea is that the late time behavior may be associ¬ 
ated wi th a geometrical effect of (relativistic) “limb dark¬ 
ening ” (jPe’eiH [20081 : IPe’er fc Rvdell20111 iLundman et al.l 
I2013T) . though a non-spheric al jet structure may be re¬ 
quired (|Deng k. Zhanal2014f) . If this is the correct in¬ 
terpretation, it implies that reliable estimation of the 
hydrodynamic parameters can be done only if the data 
are taken at the break time, or earlier. Fortunately, in 
recent years enough data are available to carry out these 


Several authors have identified a thermal emission 
component in various GRBs (detected by both the 
BATSE instrument as well as by Fermi satell ite), and 
studied their properties. As was shown by IPe’er et al.l 
(1200711 , clear identification of both the temperature and 
flux in bursts with a known redshift z (luminosity dis¬ 
tance di), can be used to infer the properties of the out- 


8 This model assumes that the acceleration is dominated by pho- 
ton pressure, and that the magnetic fields are dynamically sub¬ 
dominant. For a general treatment of t he dynamics in highl y mag¬ 
netized outflow, we refer the reader to I Gao fc Zhang] <120151) . 

9 We should emphasis the fact that a systematic study of the 
temporal behavior of GRB pulses was carried only for BATSE 
bursts. No such systematic study was carried so far for Fermi 
bursts. 
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time-dependent calculations for a substantial number of 
GRBs. 

In this work we collected the available data, and di¬ 
vided them into three categories. In category (I) we in¬ 
clude seven GRBs that fulfill the entire set of conditions: 
(1) their redshifts are known; (2) a thermal component 
was reported in the literature; and (3) a time dependent 
analysis could be performed, and therefore the values of 
r and ro were inferred. These GRBs were analyzed in 
recent years by various authors. Here, we collected the 
published data, and validated it. It case of multiple- 
pulsed GRB, we normally picked the first pulse. Two 
of these GRBs, G RB970828 and GR B990510, were de¬ 
tected by BATSE dPe’er et al.ll2007D . and the other five 
were detected by Fermi. We summarize in Table [T] the 
derived values of the parameters, as well as the references 
from which they were taken. 

In category (II) there are four GRBs detected by Fermi, 
that both: (1) show clear evidence of a significant ther¬ 
mal component, and (2) have a well defined temporal 
evolution, that enables a clear identification of the break 
time0 However, as opposed to GRBs in our category (I) 
sample, the redshifts of these GRBs are unknown. The 
dynamical properties of these GRBs can therefore be de¬ 
duced only up to the uncertainty in the redshift. The 
derived parameters of these GRBs, as given by the vari¬ 
ous authors are also shown in Table |TJ with the assumed 
values taken for the redshifts by the different authors. 
In our analysis, we did not modify the assumed redshift, 
as it was taken as the mean redshift of GRBs detected 
by the relevant instrument (note that GRB120323 is a 
short GRB). As we show below, this uncertainty does 
not affect our conclusions. 

Although the thermal component observed in the spec¬ 
tra of GRBs in our category (II) sample is statistically 
significant, it is relatively weak: in 3/4 GRBs, the flux in 
the thermal component is less than 10% of the total flux 
observed in y-rays. This poses a challen ge to the anal- 
ysis method: as wa s shown bv iZhang fe Pe'eil (|2009fl : 
lHascoet et al.1 (|2013T ) : iGao fc Zhand (120151 1. the photo- 
spheric component is suppressed if the flow is highly 
magnetized. Thus, weak thermal fluxes may be an indi¬ 
cation for highly magnetized outflow, in which case the 
dynamical calcul ations present ed in Equations [1] and [2] 
are modified fsee lGao fe ZhaneJl2015l for a full treatment 
in this case). Nonetheless, the magnetization parameter 
is unknown, and a weak thermal component does not ne¬ 
cessitate a high magnetization; it could result, e.g., from 
r p h ^ r sa t . We therefore decided to include these bursts 
in our sample. As will be shown below, the values ob¬ 
tained for both ro and T for these bursts are similar to 
those obtained for the GRBs in the rest of our sample, 
which may indicate similar dynamics. 

In category (III) we used 36 GRBs dete c ted by BATSE, 
taken from the sample of lRvde fc Pe’eil (12009I) . This is 
the largest single sample of GRBs in which a temporal 
analysis was carried out and a thermal comp onent was 
clearly identified. Out of 56 GRBs in the iRvde fc Pe’eil 
( 2009 ;) sample, we selected those GRBs in which a clear 
break time in the temporal behavior of both the tem- 

10 We omitted from our sample GRBs for which this tem¬ 
poral b ehavior could not be verified, such as, e.g., GRB100507 
IGhirlanda et abll2013l 'l. 


perature and flux was identified. Furthermore, the break 
times of the temperature and flux evolution were consis¬ 
tent with each other. Thus, although a break time was 
identified in all GRBs in our sample, the sub-sample of 
GRBs in this category is homogeneous as all GRBs are 
detected by the same instrument and identical selection 
criteria was used. 

The sample of our category (III) GRBs is given in Ta¬ 
ble [2] The values of the temperature, thermal flux and 
ratio of thermal to total flux are given at the break time; 
the derived values of T, ro and the photospheric radius, 
r p h are under the assumption of z = 1. 

While the redshifts of all the GRBs in our categories 
(II) and (III) are unknown, we point out that the addi¬ 
tional uncertainty in the estimate of ro due to the lack 
of a precise redshift is not very large. This is due to the 
fact that ro oc g?l/(1 + z) 2 , and, for the range of red¬ 
shifts typical for pre-Swift GRBs, 0.5 < z < 2.5, one 
finds 0.76 < G?L/dL( z =i) x (2/(1 + z)) 2 < 1.06. A similar 
calculation shows that the added uncertainty in the esti¬ 
mate of r is o.7 < r(*)/r (a=1) < 1 . 75 , for 0.5 < 2 < 2 . 5 . 

We estimated the additional uncertainty of GRBs in 
our category (III) due to the unknown redshifts as fol¬ 
lows E3 We used the p re-Swift d istribution of GRB 
redshifts given by l.Takobsson et all (t2006t ). We per¬ 
formed a Monte-Carlo simulation, simulating 10 6 GRBs 
drawn from this distribution, and calculated the ratios 
g?l/( 1 + z) 2 and (d^ x (1 + 2) 2 ) 1 / 4 of each GRB in our 
simulation, normalized to these values for z = 1. In 
our calculation, we assumed a standard cosmology (flat 
universe with = 0.286 and Hq = 69.6 km/s/Mpc). 
Based on the simulated results, we conclude that the av¬ 
erage values of ro and F are (ro)/ro z= i = 0.863, and 
(F)/F z= i = 1.117. The standard deviations due to the 
uncertain redshifts are cr(ro) = 0.2533 and (j(F) = 0.531. 
These values are not surprising given the above analy- 
sis and the fact t hat th e mean redshift in the sample of 
iJakobsson et al.l (12006}) is (z) = 1.345. 

The values of T and ro presented in tables nm are de¬ 
rived under the assumption of Y = 1. This is in order 
to be consistent with the results presented in the various 
references from which the data is adopted. However, in 
presenting the results in the figures below, we adopt a 
somewhat higher value of Y = 2. As explained above, 
the exact value of Y is very difficult to measure. Mea¬ 
surements based on afterglow observations reveal mixed 
results. Some works found very hi gh e ffici ency in y- 
ray production, i mplying Y < a few dCenko et alJl201ll 
iPe’er et al.l 120121) . On the other hand, several works 
found inefficient rad i ation, implying la rger values of Y 
(| Sant ana et al .1 [20~Ml : iWang et al.ll2015D . A large value 
of Y would increase the measured value of the Lorentz 
factor by tens of %, as T oc F 1//4 , while decreasing r 0 by 
a factor of a few, as ro oc Y ~ 3 / 2 . As the highest values 
of T we obtain are > 10 3 , while the lowest value of ro 
is ~ 10 6 ' 5 cm (see below), we deduce that the value of 
Y cannot be much greater than unity. In the analysis 
below, we therefore take as a fiducial value Y = 2. 

11 We omitted the full calculation for the four GRBs in our 
category (II), as this sample is not homogeneous, and the additional 
error would not affect our final conclusion. For completeness, we 
do show in Figure [T] the values of ro and T of these GRBs obtained 
for z = l. 
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TABLE 1 

List of GRBs in our category (I): known redshift, and (II): Fermi GRBs without known redshift. The derived values of r and rg in this table are done under the assumption 

Y = 1. See the text for details. 


Category Burst 

z 

T 

[keV] 

f bb 

[erg cm -2 s — 

F BB/ F tot 

r 

r 0 

[cm] 

Reference 

Comments 

I 

970828 

0.9578 

78.5 ± 4 


0.64 

305 ± 28 

2.9 ± 1.8 X 10 8 

1 

a, b 

I 

990510 

1.619 

46.5 ± 2 

7.0 X 10 -7 

0.25 

384 ± 71 

1.7 ± 1.7 X 10 8 


b 

I 

080810 

3.355 

62 

1.6 X 10 -7 

0.28 

570 ± (170) 

2.3 ± (1.2) X 10 8 

2 

c, d 

I 

090902B 

1.822 

168 

1.96 X 10 -5 

0.26 

995 ± 75 

5.2 ± 2.3 X 10 8 



I 

090926B 

1.24 

17.2 ± 1 

3 X 10 -7 

0.92 

110 ± 10 

4.3 ± 0.9 X 10 9 


e 

I 

101219B 

0.55 

19.1 ± 0.7 

8.45 ± 0.03 X 10 -8 

< 1 

138 ± 8 

2.7 ± 1.6 X 10 7 

5 


I 

110731A 

2.83 

85 ± 5 


0.47 

765 ± 200 

3.46 ± 1.1 X 10 8 

6 

d 

II 

100724B 

(i) 

38 ± 4 

2.6 X 10 -7 

0.04 

325 ± 100 

1.2 ± 0.6 X 10 7 

7 


II 

110721A 

(2) 

30 


0.09 

450 ± 200 

1 ± 0.4 X 10 7 

8 

f 

II 

110920 

(2) 

61.3 ± 0.7 


0.3 

442 ± (133) 

2 ± 1 X 10 8 



II 

120323 

(0.5) 

11.5 ± 1.5 


0.05 

145 ± (20) 

2.6 ± (0.9) X 10 9 


g 

1 IPe’er et 

al.l 120071): 2 IPaee et 

Hn 120091): 

8 IPe’er et 

al. 1 20121); ^ISerino et 

li1.1 il 20111); L. 

arsson et al. I 2015); Basha I 20131); 

7 IGuiriec 

et al.1 020111): 











Ivvani et 

al. 120131); McGlvnn &; Fermi 

GBM Collaboration! 12012); Ivvani et al. ( 20151); 

u IGuiriec et 

al. I2013j). 




a Data is based on published references. Some references omit data on Fgg, dFgg or dT. 

k GRBs 970828 and GRB990510 were detected by BATSE instruments. All other GRBs in this list were detected by Fermi-GBM. 
c Errors in drg, dT are estimated based on the data provided in the reference. 

^ In addition to Fermi-GBM, this GRB was detected by the Swift-BAT. 
e In addition to Fermi-GBM, this GRB was detected by MAXI. 

^ Errors represent uncertainty in redshift as well. 

£ Short GRB; Errors estimated from data provided in table 4 of IGuiriec et al.| H20139) . 


TABLE 2 

Sample of our category III GRBs: from IRyde fc Pe’erI < 20091) . See the text for details. 


Burst 

Trigger 

T 

[keV] 

f BB 
[erg cm -2 

s 1 ] 

F BB / F tot 

r 



r 0 

[cm] 





r ph 

[cm] 



910807 

647 

57.3 ± 2.2 

1.7 

± 

0.3 

X 

10 -6 

0.64 

257 ± 10 

6.23 

± 

1.54 

X 

10 8 

4.88 

± 

0.60 

X 

10 11 

910814 

678 

173.1 ± 12.6 

5.7 

± 

2.0 

X 

10 -6 

0.51 

549 ± 39 

8.96 

± 

4.68 

X 

10 7 

2.09 

± 

0.44 

X 

10 11 

911016 

907 

62.5 ± 2.5 

9.5 

± 

2.0 

X 

10- 7 

0.67 

247 ± 10 

4.22 

± 

1.03 

X 

10 8 

2.96 

± 

0.38 

X 

10 11 

911031 

973 

59 ± 6.9 

6.9 

± 

4.7 

X 

10 -7 

0.18 

334 ± 55 

7.00 

± 

8.33 

X 

10 7 

3.53 

± 

1.23 

X 

10 11 

920525 

1625 

148.5 ± 10.4 

1.2 

± 

0.4 

X 

10 -5 

0.43 

583 ± 42 

1.38 

± 

0.78 

X 

10 8 

4.41 

± 

0.89 

X 

10 11 

920718 

1709 

44 ± 2.3 

1.7 

± 

0.5 

X 

10 -6 

0.28 

278 ± 18 

3.02 

± 

1.66 

X 

10 8 

9.06 

± 

1.41 

X 

10 11 

921003 

1974 

15.9 ± 1.3 

7.7 

± 

3.3 

X 

10 -7 

0.59 

125 ± 12 

5.08 

± 

3.14 

X 

10 9 

2.12 

± 

0.58 

X 

10 12 

921123 

2067 

56.1 ± 1.8 

2.8 

± 

0.5 

X 

10 -6 

0.40 

306 ± 11 

4.13 

± 

1.19 

X 

10 8 

7.87 

± 

0.76 

X 

10 11 

921207 

2083 

94.8 ± 3.2 

2.0 

± 

0.1 

X 

10 -5 

0.61 

458 ± 11 

7.37 

± 

1.48 

X 

10 8 

1.11 

± 

0.09 

X 

10 12 

930112 

2127 

111.7 ± 7.5 

2.9 

± 

1.0 

X 

10 -6 

0.42 

426 ± 30 

1.13 

± 

0.63 

X 

10 8 

2.78 

± 

0.54 

X 

10 11 

930214 

2193 

100.7 ± 7.3 

6.2 

± 

2.4 

X 

10 -7 

0.81 

283 ± 23 

1.73 

± 

0.62 

X 

10 8 

1.05 

± 

0.26 

X 

10 11 

930612 

2387 

40.2 ± 4.3 

4.7 

± 

2.7 

X 

10 -7 

0.42 

208 ± 26 

4.15 

± 

3.93 

X 

10 8 

4.22 

± 

1.47 

X 

10 11 

940410 

2919 

61.4 ± 14.1 

2.3 

± 

2.7 

X 

10 -7 

0.11 

320 ± 102 

1.46 

± 

4.03 

X 

10 7 

1.94 

± 

1.26 

X 

10 11 

940708 

3067 

73.2 ± 9.8 

1.3 

± 

0.9 

X 

10 -6 

0.19 

377 ± 66 

5.51 

± 

8.25 

X 

10 7 

3.80 

± 

1.47 

X 

10 11 

941023 

3256 

46.5 ± 7.5 

3.3 

± 

1.4 

X 

10 -7 

0.61 

196 ± 29 

4.31 

± 

3.29 

X 

10 8 

2.55 

± 

0.81 

X 

10 11 

941026 

3257 

55.4 ± 6.4 

6.1 

± 

2.2 

X 

10 -7 

0.54 

233 ± 24 

3.29 

± 

2.04 

X 

10 8 

2.89 

± 

0.78 

X 

10 11 

941121 

3290 

41.1 ± 2 

1.3 

± 

0.3 

X 

10 -6 

0.34 

247 ± 14 

4.05 

± 

1.94 

X 

10 8 

8.00 

± 

1.18 

X 

10 11 

950403 

3492 

108.5 ± 4.1 

3.3 

± 

0.6 

X 

10 -5 

0.35 

598 ± 25 

3.02 

± 

1.03 

X 

10 8 

1.39 

± 

0.15 

X 

10 12 

950624 

3648 

28.9 ± 2.6 

1.3 

± 

0.7 

X 

10 -7 

0.38 

150 ± 16 

3.06 

± 

2.76 

X 

10 8 

3.06 

± 

0.89 

X 

10 11 

950701 

3658 

67.6 ± 3 

3.0 

± 

0.7 

X 

10 -6 

0.42 

334 ± 16 

3.13 

± 

1.22 

X 

10 8 

6.09 

± 

0.82 

X 

10 11 

951016 

3870 

28.2 ± 3.8 

5.9 

± 

4.7 

X 

10 -7 

0.10 

252 ± 53 

0.94 

± 

1.70 

X 

10 8 

1.17 

± 

0.48 

X 

10 12 

951102 

3891 

66.9 ± 4.9 

2.8 

± 

1.0 

X 

10 -6 

0.23 

380 ± 35 

1.28 

± 

1.00 

X 

10 8 

6.76 

± 

1.42 

X 

10 11 

951213 

3954 

52.9 ± 4.6 

9.4 

± 

4.3 

X 

10- 7 

0.22 

300 ± 34 

1.09 

± 

1.05 

X 

10 8 

4.97 

± 

1.27 

X 

10 11 

951228 

4157 

24.4 ± 0.6 

4.5 

± 

0.7 

X 

10 -7 

0.98 

128 ± 4 

3.40 

± 

0.31 

X 

10 9 

6.96 

± 

0.68 

X 

10 11 

960124 

4556 

60.6 ± 5.5 

2.3 

± 

0.9 

X 

10 -6 

0.22 

364 ± 39 

1.40 

± 

1.19 

X 

10 8 

7.29 

± 

1.71 

X 

10 11 

960530 

5478 

44.1 ± 3.6 

4.7 

± 

0.4 

X 

10 -7 

0.38 

220 ± 9 

2.55 

± 

0.43 

X 

10 8 

3.78 

± 

0.37 

X 

10 11 

960605 

5486 

64.9 ± 3.2 

1.4 

± 

0.4 

X 

10 -6 

0.54 

279 ± 15 

3.48 

± 

1.33 

X 

10 8 

3.80 

± 

0.60 

X 

10 11 

960804 

5563 

47.9 ± 3 

2.7 

± 

0.9 

X 

10 -6 

0.26 

311 ± 25 

2.93 

± 

1.98 

X 

10 8 

1.06 

± 

0.20 

X 

10 12 

960912 

5601 

68.1 ± 6.6 

7.6 

± 

3.8 

X 

10- 7 

0.33 

300 ± 34 

1.08 

± 

1.01 

X 

10 8 

2.70 

± 

0.77 

X 

10 11 

960924 

5614 

155.4 ± 8 

7.3 

± 

1.0 

X 

10- 4 

0.94 

820 ± 24 

3.17 

± 

0.48 

X 

10 9 

4.45 

± 

0.65 

X 

10 12 

961102 

5654 

62.2 ± 3.3 

8.1 

± 

2.3 

X 

10~ 7 

0.35 

283 ± 17 

1.51 

± 

0.76 

X 

10 8 

3.17 

± 

0.50 

X 

10 11 

970223 

6100 

98.3 ± 7 

3.7 

± 

1.3 

X 

10 -6 

0.32 

441 ± 35 

1.13 

± 

0.75 

X 

10 8 

4.23 

± 

0.86 

X 

10 11 

970925 

6397 

44.6 ± 2.3 

7.6 

± 

2.1 

X 

10 -7 

0.33 

243 ± 15 

2.51 

± 

1.30 

X 

10 8 

5.10 

± 

0.80 

X 

10 11 

980306 

6630 

80.6 ± 4.6 

5.5 

± 

1.1 

X 

10 -6 

0.84 

334 ± 13 

7.95 

± 

2.09 

X 

10 8 

5.83 

± 

0.82 

X 

10 11 

990102 

7293 

54.4 ± 2.2 

7.0 

± 

1.6 

X 

10 -7 

0.58 

229 ± 10 

3.85 

± 

1.18 

X 

10 8 

3.10 

± 

0.42 

X 

10 11 

990102 

7295 

100.7 ± 16.1 

1.4 

± 

0.6 

X 

10 -6 

0.68 

330 ± 45 

2.20 

± 

1.46 

X 

10 8 

1.90 

± 

0.64 

X 

10 11 


3. RESULTS 

The inferred values of the Lorentz factor T and the 
acceleration radius ro from our sample are presented in 
Figure [l] The green points represent GRBs with known 
redshift (our category I GRBs), GRBs in our category 
(II) are presented by the blue points, and category (III) 
GRBs are shown by the magenta points. The error bars 
represent statistical errors; additional errors due to the 
uncertain redsliifts are shown by the dashed (red) dots. 
The blue stars represent the values of the parameters of 
GRBs in our category (II), with assumed redshift z = 1. 
As GRBs in this category form a small fraction of our 


sample, the uncertainty in the redshift of these GRBs 
does not affect any of our conclusions. In the plot, we 
assume a fixed value of Y = 2. We further show in Figure 
[2] histograms of the values of T and r 0 . 

The most interesting result from these plots is the 
range of parameter values. We find that the range of 
values of the Lorentz factor stretches between 130 and 
1200 (with 1-cr errors inclusive), with an average value 
of (r) = 370, or (log 10 (r)) = 2.57 and standard de¬ 
viation of <7(log 10 (F)) = 0.21. These results are sim¬ 
ilar to the values of T inferred for bright GRBs by 
other method s which are in wide use, such as opac- 
ity arguments (jKrolik fc PieillTooTl iWoods fc Loeblll995l : 
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iLithwick fe Sarill200lD , the dece lerati on time of the emis¬ 
sion from__the_Jorward_shock (|Sari fe Piranl ITOOOl) : see 
also iRacusin et all (|2011l) . or the onset of the afterglow 
(iLiansr et al.ll2010l ). Furthermore, the results of Figure 
[5] indicate that the distribution of T is fairly well fitted 
by a Gaussian, suggesting that selection biases are not 
playing a major role (see further discussion below). 

The initial acceleration radius, ro, ranges between 4 x 
fO 6 < ro < 2 x fO 9 cm. The average value is (ro) = 
9.fxl0 7 cm, or (log 10 (r 0 )) = 7.96 and standard deviation 
of cr(log 10 (r 0 )) = 0.56. Similarly to the distribution of T, 
the distribution of log 10 (ro) is also close to a Gaussian 
(see FigurcOQ The average value found corresponds to 
minimum variability time of St = ro/c ~ 3 x 10 -3 s, and 
is about 30 times the gravitational radius of 10 M 0 black 
hole. While this value is fully consistent with our data, as 
all the GRBs in our sample are long, and their variability 
time is much longer than 10 ms, this value is different 
from the commonly assumed value of the acc ele ratio n 
radius, » few gravita tional radii (e.g., fMcszaros! 120061 
iKumar fc Zhandl20l4 and references therein). 

A priori, we do not expect any correlation between 
the values of T and of ro, as these are two independent 
parameters of the fireball model. Nonetheless, it is inter¬ 
esting to search for such a correlation, since if it exists 
it could put interesting constraints on the nature of the 
progenitor. This is due to the fact that if indeed ro is 
related to the gravitational radius, than ro oc M, where 
M is the mass of the central object, while T oc L/Me 2 , 
where M is the mass ejection rate. 

Linear fitting (with points weighted by the errors in 
both directions) reveals best fitted values of log 10 (r) = 
ai + ci 2 log 10 (ro), with a\ = 3.31 ± 0.66 and «2 = 
—0.10 ± 0.08. The correlation is thus very weak, with 
intrinsic scatter 0.44 ± 0.05 (Spearman rank -0.31 with 
significance of 0.035, equivalent to 2.1er). While this 
analysis is done on a non-homogeneous sample, a sim¬ 
ilar analysis carried out on GRBs in our category (III) 
only, which do form a homogeneous sample, is not sig¬ 
nificantly different. We can conclude that our data does 
not points to any correlation between the values of ro 
and T. 

4. DISCUSSION 
4.1. Possible selection bias ? 

The derived values of T span about one order of magni¬ 
tude, between 10 2 —10 3 . This range is similar to the typ¬ 
ical values of T inferred from other methods, as discussed 
above. The values of ro span a wider range of roughly 
two orders of magnitude, between 10 65 —10 9 ' 5 cm. These 
values are higher than the typically assumed value of 
ro < 10' cm. 

The values of T and ro are derived from the iden¬ 
tified values of the observed temperature, T ob ■ as well 
as the fluxes Ff b fj and F° b t . Typical values of the ob¬ 
served GRB fluxes a re in the range of F° b f ; ~ 10~ 7 — 
10~ 4 erg cm~ 2 s -1 (iKaneko et~ahl 120061 iGruber et al.l 
12014 Ivon Kienlin et al.ll20l4 see also Tables ED©. The 
ratio of thermal to total flux cannot be less than a few % 
(typically, it is of the order of 50 %, see Table [2]), other- 

12 We note that as about 75% of the GRBs in our sample are in 
category (III), they have the dominant effect on the distributions. 



Fig. 1.— Lorentz factor, T, vs acceleration radius, ro, of GRBs 
in our sample. Green points are GRBs in our category (I); Blue 
points are GRBs in our category (II), while magenta points are 
GRBs from our category (III) (see text for details). Solid error 
bars represent statistical errors, while dashed error bars represent 
additional uncertainty due to unknown redshifts of GRBs in our 
category (III). The stars show the location of the parameters of 
GRBs in our category (II) with assumed redshift z = 1. A linear 
fit reveals T a rj with a = —0.10 ± 0.09 and a very weak cor¬ 
relation. The values of T and ro presented are derived under the 
assumption Y = 2. For a different value of Y, the results can be 
scaled according to T oc F' 1 / 4 and ro oc Y ~ 3 / 2 . 



Fig. 2.— Histograms of the mean values of log 10 (ro) (left) and 
logio(r) (right). Blue are for the enti re samp le , while mag enta are 
for GRBs in our category (III) (the IR.vde fc Pe’eri 120091 1 sample 
only, and green are for GRBs in our categories (I) and (II). 

wise the thermal component could not have been clearly 
identified. Furthermore, the BATSE detector is sensitive 
at the range 20 keV - 2 MeV, while the GBM detector has 
even broader sensitivity range, 8 keV - 40 MeV. Thus, a 
thermal peak could be identified for the BATSE bursts 
at temperatures 10 3 < T ob < 10 6 eV. Using these obser¬ 
vational constraints, we plot in Figure [3] the range of ro 
and r that could have been derived from observations of 
temperature and fluxes in this range using Equations [Q 
O under the assumption of z = 1. 

Clearly, the possible range of values of both T and ro 
are much greater than the observed ones. Values of T are 
possible in the range 10 15 < F < 10 3 5 , and those of ro 
could span a much broader range, 10 5 5 < ro < 10 13 cm 
(see Figurc[3]). The addition of uncertainty in the redshift 
would imply that, in principle, an even broader range 
of parameter values could be obtained. However, this 
is not observed: the inferred values of T and ro span 
only a small part of the range allowed by the detector’s 
capabilities. Combined with the fact that the histograms 
of both F and r 0 are close to Gaussian, we thus conclude 
that selection effects are likely to be ruled out. 

Our sample is not homogeneous, as we use data ob¬ 
tained by two different instruments (Fermi and BATSE 
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T [eV] 

Fig. 3. — Derived values of the Lorentz factor, T (solid), 
and ro (dashed lines), as a function of the observed tem¬ 
perature and flux within typical detector’s capabilities. The 
lines, from top to bottom, represent different GRB fluxes: 
10 -4 ,10 -5 ,10 -6 ,10 -7 erg cm -2 s -1 (blue, greed, red, magenta, 
respectively). Lines derived for Fbb / Ftot = 1/2 and z = 1. The 
stars show the location of the extremum observed values of T, and 
the squares the extremum observed values of ro- 

detector). Furthermore, for the GRBs in our category 
(I) and (II) we rely on analysis carried by various au¬ 
thors. Thus, there is a risk of bias in the results. This 
risk does not exist for GRBs in our category (III), as all 
the data is taken from the BATSE instruments, and all 
the analysis was carried by us. Our category (III) GRBs 
constitute the largest fraction of GRBs in our sample 
(36/47 ~ 75%). When conducting a separate analysis to 
the data in our category (III) GRBs only, we obtain sim¬ 
ilar results to those obtained when analyzing the data in 
categories (I) and (II) (see Figures [T| and [2j, indicating 
that selection bias, if any, do not significantly affect the 
results. 

An additional source of uncertainty is the unknown 
value of Y > 1, the ratio between energy released in 
the explosion and the energy observed in 7-rays. As 
discussed above, the value we chose here, Y = 2, is a 
realistic approximation, based on (1) estimates from the 
literature; (2) the fact that the values of T obtained us¬ 
ing this value are consistent with measurements based 
on other methods; and (3) the fact that a higher value 
would imply, in some cases, values of ro lower than the 
gravitational radius, which is unphysical. Nonetheless, 
we stress that as the derived values of both T and ro de¬ 
pend on Y they are obviously sensitive to the uncertainty 
in the value of this parameter. 

4.2. How ubiquitous is thermal emission in GRBs ? 

The size of the sample considered in this work, 47 
GRBs, is a small fraction of the total number of GRBs 
observed to date. All GRBs in our sample are character¬ 
ized by a clearly detected, significant thermal (Planck- 
like) component, which shows a well defined, repetitive 
temporal behavior. We point out that a necessary condi¬ 
tion for performing this analysis is the ability to carry out 
a time-resolved spectral analysis. This limits the num¬ 
ber of GRBs in which this analysis could be executed 
to only those that show relatively long, smooth pulses, 
and, in addition, are bright. Furthermore, we point out 
that the traditional use of a “Band” function on its own 
in fitting the data, excludes the possibility of identifying 
a thermal component on top of a non-thermal one, as 


the “Band” function is simply not capable of capturing a 
thermal peak, where it exists. Thus, in order to obtain a 
reliable estimate of the fraction of GRBs with a thermal 
component, a re-analysis of the en t ire sample is needed . 
Recent works (jBurgess et al.l[20H : iGuiriec et, al.ll2015bl f 
in which re-analysis of the data was done, indicate that 
tens of % of bright GRBs may show evidence for the ex¬ 
istence of a distinct thermal component. As shown in 
Table ED GRBs which do show evidence for a thermal 
component span a wide range of properties in terms of 
fluxes, redshifts and light curves, further supporting the 
idea that they are ubiquitous. 

The existence of a thermal component can be fur¬ 
ther deduced indirectly, in those GRBs in which a clear 
“Planck” spectrum is not observed. This is done by ana¬ 
lyzing the low energy spectral slopes obtaine d usin g the 
“Band” fits. As was shown already in 1998 (IPreece et all 
119981120021 IGhirlanda et all [200311 . the low energy spec¬ 
tral slopes of over 85% of the GRBs are too hard to 
be consistent with the (optically thin) syn chr otro n emis¬ 
sion. A more recen t analysis (lAxelsson fc Borgonovol 
120151 lYu et al]l2015bj) reveals that the spectral width of 
nearly 100% of GRB pulses are inconsistent with having 
a synchrotron origin, unless additional assumptions are 
made, such as variation of the magnetic field with radius 
dPe’er fe Zha.ngl [20061 lUhm fc Zhang|l2014t iZhang et al.1 
1201511 . 

One possible solution that can explain the observed 
spectra of those GRBs that do not show clear evidence of 
a distinct thermal component (which are the vast major¬ 
ity) is that they do in fact originate from a thermal com¬ 


ponent, tha t is modified due to sub-photospheric energy 
dissipa tion (fPe ’er et al.l[2005l [20061: lGlanniosll2008L ( 2012: 


Vurm etaTI 20 lT Lazzati et, all 20131 Deng fc Zhand 


^Ol^jKere^^Levinson 2014 iGao fc Zhand 
Chhotrav fc Lazzatil 120151 IVu rm fc Beloborodov! 


Ahlgren et al.112015T ). If this hypothesis is correct, then 


the origin of the “Band” peak is Comptonization of 
the thermal component. A similar mechanism may 
be responsible for the non-thermal part of the spectra 
in those GRBs in which a thermal component is pro¬ 
nounced, though an alternative synchrotron model for 
the non-thermal part is also a possibility in these bursts 
dBurgess et alJ20fHl20l4lYu et al.ll2015alf . Thus, in this 
model, the key difference between GRBs in our sample 
(in which a thermal component is identified) and other 
GRBs, is that for GRBs in our sample the “Planck” 
part of the spectrum suffers only relatively weak dis¬ 
tortion. This could be due to the lack of significant 
dissipation at small (sub-photospheric) radii. We point 
out that in the classical fireball model, where internal 
shocks may be the main source of energy dissipation, 
there is only a low limit constraint on the location of 
the internal shocks, r > V 2 cSt , where St is the vari¬ 
ability time. Thus, for smooth bursts, for which St is 
relatively long, i ntern al sho cks o ccur far above the pho¬ 
tosphere dRees fc Meszarosll2005H . As a result, if indeed 
internal shocks are the main cause of kinetic energy dis¬ 
sipation, one naturally expects a pure photospheric com¬ 
ponent to be more pronounced in smooth bursts. We 
stress, though, that so far no systematic study of pos¬ 
sible correlation between the smoothness of lightcurves 
and the existence and properties of a thermal component 
was carried. 
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An alternative possibility of explaining the var¬ 
ious s pectra is due to diffe rence in ma g netiza¬ 
tion (fZhang fc Pe’eH 120091 : iZhang fe Yanl 2011 : 
Beniamini fc Piranl 120141 : iGao fc Zhang! 2015; 

Begue fe Pe erl I2015T) . In highly magnetized out¬ 
flows (characterized by er = ^ 1, where Em is 

the energy stored in the magnetic field and Ek is the ki¬ 
netic energy), photospheric emission is both suppressed 
(e.g., IZhang fe Pe’eH I2009D and is further shifted to 
higher energies due to suppression of photon production 
below the photosphere. Thus, according to this model, 
the difference between the GRBs in our sample and 
other GRBs is that the outflow of the GRBs in our 
sample happen to have a weaker magnetization, a < 1. 
However, Currently, there is no clear way to determine 
the magnetization of GRB outflows without additional 
assumptions. On the other hand, we point out that an 
additional assumptions, e.g., about the value of ro can 
be used to interpret the data in the framework of the 
magnetized outflow model, and use this interpretation 
to infer the value_of the magnetization parameter (see 
IGao fe ZhangH2015l for details). 

Finally, we point out that the larger part of our sample 
(our category III) are GRBs detected by BATSE. This 
is mainly due to the fact that a systematic analysis in 
search for a thermal component was carried out, so far, 
only on BATSE data, but not on data obtained by Fermi. 
It is therefore not clear yet what fraction of GRBs de¬ 
tected by Fermi show clear evidence for th e existence of 
a ther mal component. A recent analysis (SGuiriec et all 
l2015al) indicates that the fraction of thermal compo¬ 
nent in Fermi bursts may be similar to those observed 
in BATSE bursts, though further analysis is still needed. 

4.3. Possible evidence of a massive progenitor 

Wi thin the framework of the “collapsar” model 
(lWooslevlll993D . the jet drills its way through the collaps¬ 
ing star. For a freely expanding outflow, it is typically 
assumed that the initial acceleration radius, ro is close to 
the gravitational radius of the newly-formed black hole, 
namely ro < 10' cm. However, the results presented 
here indicate a rather different value of (ro) = 10 s cm, 
namely, about a factor 10 higher. 

We should stress that the idea that r 0 should be « few 
gravitational radii is of theoretical origin, and is correct 
in the classical (low magnetization) “fireball” model, in 
the limit of free outflow. A small acceleration radius is 
associated, via light crossing time argument, with short 
variability timescales, as 5t = r^/c. However, the min¬ 
imum observed variability time in GRB s is > 10 ms, 
with an average value of « 500 ms (Golkho u fe Butlen 
I2014D . Variability timescales of 10 ms correspond to 
ro ~ 10 85 cm. Thus, in fact, there is no observational 
evidence supporting values of ro much smaller than this 
value. 

Adopting the basic framework of the “collapsar” 
model, the collimated outflow within the collapsing star 
cannot be described as a free outflow. As the jet prop¬ 
agates through the star, it forms a hot “cocoon” made 
of hot stellar material heated by the jet itself that sur¬ 
rounds the jet. Furthermore, as the jet drills a funnel 
through the stellar material, it is confined by the fun¬ 
nel walls. The jet thermal pressure decreases with dis¬ 
tance from the origin. Thus, if the external pressure 


decreases slower than the jet pressure, a recollimation 
shock must form. The shape of this recollimation shock 
depends on the balance between the external pressure 
and the iet ram pressure, a nd can be c alcul ated ana¬ 
l ytica lly (IKomissarov fe Falli 119971 : INalewaiko fe Sikoral 
f2009l) . 

This recollimation shock is very efficient i n dissipa ting 
the kinetic energy (INalewaiko fe Sikorall2009{ : INalewaiko! 
I2012D . Thus, while the Lorentz factor increases up to the 
recollimation shock radius, there is a sharp drop in the 
outflow velocity as it encounters the recollimation shock 
to r > 1, before the outflow re-accelerates above this 
radius. This is clearly seen in various numer ical models 
of jet propagation inside stellar cores (e.g.. lAlov et al.1 
l2002t IMorsonv et al.ll2007t IMizuta fe Alovll2 009) as well 
as in analytical modeling (IBromberg et alJl201lD . The 
radii of the recoll imation shocks are typically seen at 
~ few x 10 9 cm (IMizuta fe Alovl I2009t Mizuta fe Iokal 
120131 iLopez-Camara et al.ll2013D though in some of the 
simulations, recollimation shocks can be identified only 
above 10 9 cm due to lack of numerical resolution at 
smaller radii. 

Once the jet completes its crossing of the stellar core, 
both the jet and the cocoon expand into the stellar en¬ 
velope and into the interstellar medium (ISM). At this 
stage, the external pressure rapidly drops, and the ra¬ 
dius of the recollimation shock gradually increases, until 
the shocks eventually disappear and the flow becomes 
free. However, this stage lasts typically for a dura¬ 
tion of > te n seconds (at least several sound crossing 
times); (e.g.. IMorsonv et al.ll2007t IMizuta fe Alovl 120091 : 
ILopez-Camara et al.ll2013l~ During this epoch, the rec- 
ollimation shock is roughly at its initial location. It can 
be shown that the presence of a recollimation shock does 
change the length scale of the jet (ro) in certain con¬ 
ditions, especially if the jet is brought back to causal 
contact by the collimation. In that case, the measured 
value of ro would not correspond to the size of the accel¬ 
eration region nor to the radius at which the shocks hap¬ 
pen, but to roTsh, where r s h is the Lorentz factor of the 
shocked jet. In our case, a Lorenz factor of ~ 30 would 
be able to reconcile the observations with the theoreti¬ 
cal prediction of a small acceleration region of ~ 10 7 cm. 
While simulations show that typically r s h < 30 (see, e.g., 
IMorsonv et al.l[2007t iLopez-Cdmara et alj|2013l ) , they do 
also show that shocks cause mass entrainment on the jet 
and that the jet pro pagati o n is no n-co nical at least for a 
few tens of seconds ([Morsonv et alJl2007f ) . Both these ef¬ 
fect contribute also to causing a measured value of ro that 
is larger than the size of the acceleration region (Lazzati, 
Pe’er, Ryde, et al in preparation). 

We may thus conclude that the results obtained here, 
of r 0 ~ 10 7 — 10 9 cm (despite all the uncertainties 
discussed above) are consistent with the picture that 
emerges from the numerical simulations, and that they 
provide indirect evidence that even the brightest long du¬ 
ration bursts at cosmological redshifts can be produced 
by the core-collapse of massive stars. Furthermore, this 
implies that if indeed the progenitors of short GRBs are 
merger of binaries rather than collapses of massive stars, 
then the values of r o are expected to be much smaller in 
these bursts. While one of the bursts with a high value of 
ro in our sample (GRB120323) is categorized as short, we 
note that the discrimination between “long” and “short” 
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is statistical in nature, and thus it is possible that this 
burst may originate from a collap se o f a massive st ar, 
similar to classical long GRBs Isee lZhang et al.ll2009D . 

5. SUMMARY AND CONCLUSIONS 

In this work, we analyzed a sample of 47 GRBs which 
show significant thermal emission component. We split 
our sample into three categories: (I) 7 GRBs with known 
redsliifts; (II) 4 GRBs detected by Fermi, in which the 
redshifts are unknown; and (III) 3 6 GRBs detec t ed by 
BATSE, selected from the sample of IRvde fc Pe’eil (120091 ) 
based on clear identification of the break time. We ana¬ 
lyzed the dynamical properties of these GRBs based on 
standard “fireball” mo del as sumptions using the method 
derived inlPe ’er et al.l (120071) . We found a mean Lorentz 
factor of (r) = 370, consistent with values obtained by 
different methods. However, we further found that the 
acceleration radius ro ranges between 4 x 10 6 < r o < 
2 x 10 9 cm, with mean value (ro) — 10 8 cm. This is 
~one orders of magnitude above the commonly assumed 
value of ~ 10' cm. We found only a very weak corre¬ 
lation between the values of T and r o, that cannot be 
used to put further constraints on the properties of the 
progenitor. We showed that the derived values of ro are 
only weakly sensitive to the uncertainty in the redshifts. 

We showed in 44.1l that selection biases are unlikely to 
affect our obtained results. We further argued f fl4.2l) that 
despite the fact that GRBs with a pure thermal emission 
component are relatively rare, the existence of such a 
component is likely to be very ubiquitous, and possibly 
it exists in nearly all GRBs, albeit in most of them it 
is distorted. We claimed that within the context of the 
standard fireball model, it is more likely to detect a pure 
thermal component in long, smooth GRBs, that are less 
variable, and as a result internal shocks occur above the 
photosphere. Finally, we argued in t)4.3l that the most 
likely interpretation of the values of ro are propagation 
effects in massive progenitor stars, such as collimation 
shocks, mass entrainment, or non-conical expansion. 

An alternative scenario that can explain the relatively 
large values of ro is that the acceleration in fact begins 
at radii much greater than the gravitational radius of the 
newly formed black hole. This could be due to initially 
large jet opening angle, followed by a strong poloidal 
recollimation at ~ 10 2 — 10 3 r g that accelerates the jet. 
Some evidence for a similar sce nari o exists i n a d ifferent 
object, the active galaxy M87 (IJunor et al.lll999h . This 
scenario requires very strong poloidal magnetic fields, 
whose existence in a GRB environment is uncertain. 

The calculation of r 0 and T in Equations [Tj [2] were 
done under the assumption of the “classical” (weakly 
magnetized) “fireball” model. The values of T that we 


found are consistent with the values found using other 
methods, that are independent of some of the underlying 
“fireball” m odel assu mptio ns, su ch as afterglow measure¬ 
ments fe.g.. lRacusin et aHl201lD . These values are more 
than an order of magnitude higher than the values of T 
inferred in active galactic nuclei (AGNs) and X-ray bi¬ 
naries Th us, w hile magn etic acceleration (e.g., via the 
iBlandford fe Znaieklll977l process) is likely to occur in 
XRBs and AGNs, currently there is no evidence that a 
similar mechanism is at work in GRBs as well. 

In addition to the GRBs considered in this work, 
thermal emission was found recently in several GRBs 
in the X-ray band, extending well into the afterglow 
phase. These i nclude both low luminosity GRBs such 
as GR B060218 (|Campana et al .1120061 iShcherbakov et al.1 
I2013D or GRB100316D (jStarling et alJ 120111) but, also 
many GRBs with typ i cal lu m inosities (e.g.. [Page et al.1 
201 it iStarling et al. I2012t ISoarre fc Starling! 2012t 
Friis fe Watsonl l2013t iBellm et al.l 12014 iSchulze et al.l 


2014 iPiro et al. I2014D . Common to all GRBs in this 


category are (1) the fact that the thermal component is 
observed well into the afterglow phase; and (2) the in¬ 
ferred values of the Lorentz factors are at least an order 
of magnitude lower than that of GRBs in our sample: 
r < few tens, and in some cases much lower, T > 1. 
These results indicate that most likely the physical ori¬ 
gin of the thermal component in bursts in this sample is 
different than bursts in our sample, the leading models 
being supernovae shock breakout and emission from the 
emerging cocoon. Given this different origin, we excluded 
these bursts from the analysis carried out here. 

Finally, we stress that as the calculations here are car¬ 
ried under the standard assumptions of the classical “fire¬ 
ball” model, any inconsistency that may be found in the 
future between the values of tq derived by our method 
and those derived in alternative ways (e.g., via variability 
time argument) would question the validity of the entire 
“fireball” model. 
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